MULTIWAVELET PACKETS AND FRAME PACKETS OF ^^(R'^) 



BISWARANJAN BEHERA 

Abstract. The orthonomial basis generated by a wavelet of L^(M) has poor fre- 
quency locahzation. To overcome this disadvantage Coifman, Meyer, and Wicker- 
hauser constructed wavelet packets. We extend this concept to the higher dimen- 
sions where we consider arbitrary dilation matrices. The resulting basis of L^(M'^) 
is called the multiwavelet packet basis. The concept of wavelet frame packet is also 
generalized to this setting. Further, we show how to construct various orthonormal 
bases of L'^{R'^) from the multiwavelet packets. 



1. Introduction 

Consider an orthonormal wavelet of L^(]R). At the j-th resolution level, the 
orthonormal basis {ipjk '■ j,k & Z} generated by the wavelet has a frequency local- 
ization proportional to 2^. For example, if the wavelet ip is band-limited (i.e., ^ 
is compactly supported), then the measure of the support of ijpjk)'^ is 2-' times the 
measure of the support of since 

where 

V-ifc = 2^^/V(2^' ■ -fc), j,fcez. 

So when j is large, the wavelet bases have poor frequency localization. Better 
frequency localization can be achieved by a suitable construction starting from an 
MRA wavelet basis. 

Let {Vj- : j G Z} be an MRA of L^(M) with corresponding scaling function and 
wavelet ip. Let Wj be the corresponding wavelet subspaces: Wj = splipjk : G Z}. 

In the construction of a wavelet from an MRA, essentially the space Vi was split 
into two orthogonal components Vq and Wq. Note that Vi is the closure of the linear 
span of the functions {25(^(2 • —k) : k G Z}, whereas Vq and Wq are respectively the 
closure of the span of {</?(■— /c) : k} and {■?/'(■ — /c) : k}. Since ip{2-—k) = Lp (2(- — |)), 
we see that the above procedure splits the half-integer translates of a function into 
integer translates of two functions. 

In fact, the splitting is not confined to Vi alone: we can choose to split Wj, which 
is the span of {ip{2^ ■ —k) : k} = {ip {2\- — ^)) : fc}, to get two functions whose 
2-0'~i)/c translates will span the same space Wj. Repeating the splitting procedure 
j times, we get 2^ functions whose integer translates alone span the space Wj. If 

Key words and phrases, wavelet, wavelet packets, frame packets, dilation matrix. 
Supported by a grant from The National Board for Higher Mathematics, Govt, of India. 

1 



2 



BISWARANJAN BEHERA 



we apply this to each Wj, then the resulting basis of L^(]R), which will consist of 
integer translates of a countable number of functions (instead of all dilations and 
translations of the wavelet will give us a better frequency localization. This 
basis is called "wavelet packet basis". 

The concept of wavelet packet was introduced by Coifman, Meyer and Wicker- 
hauser ||^, ^. For a nice exposition of wavelet packets of L^(]R) with dilation 2, see 
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The concept of wavelet packet was subsequently generalized to by taking ten- 
sor products 0. The non-tensor product version is due to Shen |T^. Other notable 



generalizations are the biorthogonal wavelet packets ||^, non-orthogonal version of 
wavelet packets the wavelet frame packets on M for dilation 2, and the orthog- 
onal, biorthogonal and frame packets on M'^ by Long and Chen [|l^ for the dyadic 
dilation. 

In this article we generalize these concepts to for arbitrary dilation matrices 
and we will not restrict ourselves to one scaling function: we consider the case of 
those MRAs for which the central space is generated by several scaling functions. 

Definition 1.1. A d x d matrix A is said to be a dilation matrix for M*^ if 

(i) C Z'^ and 

(ii) all eigenvalues A of A satisfy |A| > 1. 

Property (i) implies that A has integer entries and hence |dety4| is an integer, 
and (ii) says that | det A\ is greater than 1. Let B = A^, the transpose of A and 
a = I dety4| = | deti?|. Considering Z*^ as an additive group, we see that AZ''- is a 
normal subgroup of Z"^. So we can form the cosets of AZ'^ in Z"^. It is a well known 
fact that the number of distinct cosets of AZ'^ in Z'^ is equal to a = | dety4| ([^ 



A subset of Z which consists of exactly one element from each of the a cosets of 
AZ"^ in Z'^ will be called a set of digits for the dilation matrix A. Therefore, if Ka 
is a set of digits for A, then we can write 

Z''= [j (AZ^' + fi), 

where {AZ''- + /i : G Ka} are pairwise disjoint. A set of digits for A need not be 
a set of digits for its transpose. For example, for the dilation matrix M = ( ? g) 
M^, the set { (g) , (q) } is a set of digits for M but not for M*. It is easy to see that if 
K is a set of digits for A, then so is K — fi, where fi ^ K. Therefore, we can assume, 
without loss of generality, that E K. 

The notion of a multiresolution analysis can be extended to L^(M'^) by replacing 
the dyadic dilation by a dilation matrix and allowing the resolution spaces to be 
spanned by more than one scaling function. 

Definition 1.2. A sequence {Vj : j E Z} of closed subspaces of L^(M'^) will be 
called a multiresolution analysis (MRA) of L^(M'^) of multiplicity L associated with 
the dilation matrix A if the following conditions are satisfied: 
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(Ml) Vj C Vj+i for all j G Z 

(M2) Uj^zVj is dense in L'^{W^) and DjezVj = {0} 

(M3) / e Vj if and only if f{A-) e Vj+i 

(M4) there exist L functions {y^i, <^2, • • • , V^l} in Vq, called the scaling functions, 
such that the system of functions {ifi{- — k) : 1 < / < L, G Z'^} forms an 
orthonormal basis for Vq. 

The concept of multiplicity was introduced by Herve |]12| in his Ph.D. thesis. 
Since {^Pi{- — k) : 1 < I < L, k E Z"'} is an orthonormal basis of Vq, it follows from 
property (M3) that {a^^VK^^' ■ -k) : 1 < / < L, G Z^} is an orthonormal basis 
of Vj. Observe that if / G ^^(M^), then 

{a^^^fiA^ ■ -k)yiO = a-^/^e-'^^~'^^^) f{B'^i), i eR\ ke Z^. 
The Fourier transform of a function / G -/^^(M'^) is defined by 



To define the Fourier transform for functions of L^(R'^), the operator T is extended 
from n L^(R°'), which is dense in L'^{W^) in the L^-norm, to the whole of L^(R'^). 
For this definition of the Fourier transform, Plancherel theorem takes the form 

First of all we will prove a lemma, the splitting lemma (see which is essential 
for the construction of wavelet packets. We need the following facts for the proof of 
the splitting lemma. 

(a) Let T'^ = [-vr, vr]'^ and / G L'^{R'^). Since R'^ = U^ez^^'^ + 2A;7r), we can write 



/ { E fix + 2k7r)\dx. 



(1) / fix)dx 

(b) Let {sk '■ k G Z''} G /^(Z'^) and Kb be a set of digits for the dilation matrix 
B. As Z'^ can be decomposed as Z'^ = U^g/^g(5Z'^ + /i), we can write 



(2) E = E E ^M+Bfc- 

(c) Let Kb be a set of digits for B. Define 

Qo= U 5-i(T'^ + 2/i7r). 

Since i^^s is a set of digits for B, the set Qo satisfies y^kei,<^{.QQ + 2A;7r) = R'^. This 
fact, together with IQol = (2vr)'^, implies that {Qq + 2A;7r : A; G Z'^} is a pairwise 
disjoint collection (see Lemma 1 of |10|). Therefore, 

(3) f f{x)dx= f I ^ /(x + 2A;7r)|(ia;, for / G ^^(R''). 
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A function / is said to be 27rZ'^-periodic if f{x + 2kn) = f{x) for all k E Z'^ and 
for a.e. x E R.'^. 

2. The Splitting Lemma 

Let {yDi : 1 < I < L} he functions in L'^{R'^) such that {<^i{- - k) : 1 < I < L, k e 
U^} is an orthonormal system. Let V = 'Wp{a^^'^(pi{A ■ —k) : /, k}. For 1 < I, j < L 
and < r < a — 1, suppose that there exist sequences {hlj/^ : k G Z'^} G P{Z'^). 
Define 

(4) f[{x)=j: Eh^iy/'^fM^-k). 



Taking Fourier transform of both sides 

i=i fegZ'* 

(5) = th^B-'OvAB-'O, 
where 

(6) h\.iO = E a-i/^/i[^.,e-«'^), I < I, j < L, < r < a - 1, 

and h^j is 27rZ'^-periodic and is in L^(T'^). Now, for < r < a — 1, define the L x L 
matrices 



(7) /7,(0 = (/^I^(O)^^ 
By denoting 

(8) $(x) = 

(9) $(0 = (^i(0,---,^l(0)*, 

we can write (|) as where F,(a;) = (/[(x), /^^(x), . . . , and ^,(0 = (/[(O, (0, • • • , 

The following well known lemma characterizes the orthonormality of the system 
— k) : 1 < I < L, k E Z'^}. We give a proof for the sake of completeness. 

Lemma 2.1. The system — k) : 1 < I < L, k E Z"^} is orthonormal if and 

only if 



Y: + 2k7r)0i{^ + 2k7T) = 6ju 1<3,1<L. 

Proof: Suppose that the system {ipi{- — k) : 1 < I < L, k E Z'^} is orthonormal. 
Note that {ipj{- - p),(pi{- - q)) = {iPj,(pi{- - (q-p))) for 1 <j,l <L and p,qe Z'^. 
Now 
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Therefore, the 27rZ'^-periodic function Gji{^) = Yl + 2/c7r)(^z(^ + 2A;7r) has 

Fourier coefficients Gji{—p) = Sji6op, p E Z'^ which imphes that Gji = 6ji a.e. By- 
reversing the above steps we can prove the converse. □ 

Let M*(^) be the conjugate transpose of the matrix M(^) and 1^ denote the 
identity matrix of order L. 

Lemma 2.2. (The splitting lemma) Let {^pi : I < I < L} be functions in L^(M'^) 
such that the system {a}/'^Lpj{A ■ —k) '■ 1 < j < L, k E 7/} is orthonormal. Let V 
he its closed linear span. Let K be a set of digits for B. Also let f[, Hr be as above. 
Then 

{f[{- -k): < r < a - 1, 1 < I < L, k E Z'^} 
is an orthonormal system if and only if 

(10) E Hr{^ + '^B-'i^tt)H*{^ + 25-V7r) = < r, s < a - L 

Moreover, {f[{- ~ k) : < r < a - 1, 1 < / < L, kEZ'^} is an orthonormal basis 
of V whenever it is orthonormal. 

Proof: For 1 < /, j < < r, s < a — 1 and p E Z*^, we have 

(/;,//(■ -p)> 



1 



(2vr) 
1 

W) 



(27r) 



[ EE h^,rniB-'OhL{B''O0rn{B-'OVniB-'Oe'(P'^^dC 

(by i)) 

Kz [ E EE {h^,m{B-\^ + 2k7r))hUB~\^ + 2kn)) 



■^UB-\^ + 2kn))if^iB-\^ + 2fc7r))| e^<P.«+2Mrf^ (by (g) 
JT^J EEE h^B-'^ + 2i?-Vvr)/^?Ji?-^e + 25-Vvr) 



^ ^^{B-\i + 2/i7r) + 2kTi)^n{B-\i + 2^nx) + 2kTx) 
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(by @) 

•5„„e^<P'«>rf^ (by Lemma ^) 

Therefore, 

L 



^ E E f^jmiB-'^ + 25- V7r)/iL(i?-^e + 25-Vvr) = for a.e. ^ G 

/xe_fsr m=i 

L 



nd 



^ E E^-^(e + 25-V^)/iL(e + 25-V7r) = for a.e. e e 

^ Hr{i + 2B~^fnT)H*{^ + 25-V7r) = 5r Jl for a.e. ^ G R'^. 

We have proved the first part of the lemma. Now assume that {fi{- — k) : 
< r < a - 1, 1 < / < L, k eZ'^} is an orthonormal system. We want to show 
that this is an orthonormal basis of V. Let f E V. So there exists {cjp : p G Z^} G 
l<j<L such that 

L 

/(a;) = E E Cj>a^/Vj(^a; -p). 

Assume that / ± — k) for all r, /, k. 
Claim: / = 0. 

For all r, Z, /c such that 0<r<a — 1, 1 < I < L, k E Z'^, we have 
= - A:), /) = {f[{- - fc), E E c,paVV,(A ■ -p) 

-1/2 PL L , _^ , 

= 7^/ E/^Lli^-^O^^l^-^Oe-^^'^'^^E E9?e^^^"^'"^^.(i?-^e)rfe 

l^TTJ JlRd m=l j=lpeZ'* 

(by ®) 



(27r)' 



/ E hUOMO E E c-^,(Oe-^<'='^^>e^<^'«>rfe U - 

JlRd m=l i=lpeZ<* 
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7^ / EE h\^{i + 2gvr)^„(e + 2g7r) 



■ E E cJ^vA^ + 2g7r)e-^<'='^(«+2''-)>e^<P'«+2''-)de (by (g)) 



7^ / E E E hUOcfpi E 0m{^ + 2gvr)^,(e + 2g7r)| 

7^ / E E /^L(Oc;^e-<'='^«)e^<^'«)rfe (by Lemma O 

7^ E / EE /^L(Ocwe-^<'='^^^e^<^'«)rfe 

7^ E / E E /^L(e + 25-Vvr)c^e-K'='^(«+^^-V.)) 
/ (e E E/^L(e + 2S-Vvr)c^e^<^'«+^^"^-)| 



(2vr) 



Since | (2-i)<i ^~^^'''^'^ '■ k Z'^^ is an orthonormal basis for L'^{B~^T'^), the above 
equations give 

E E E c^e*^«+'''"'^"'^)/iL(^ + 25-Vvr) = a.e. for all r, /. 



For m = 1, 2, . . . L, define 

(11) C^„^(0 = E c^pe-<€'f>. 

So we have 



(12) ^ E C„(^ + 25-V7r)/iL(e + 2S-Vvr) = 0, 0<r<a-l, 1</<L. 

/xGA' m=l 

Equations ([T0|) are equivalent to saying that for < r < a — 1, 1 < I < L and 
for a.e. ^ G M'^, the vectors 



(^/i[^(e + 25-Vvr) : 1 < m < L, ^ G 



are mutually orthogonal and each has norm 1, considered as a vector in the aL- 
dimensional space C'^'^, so that they form an orthonormal basis for C""^. Equation 
(P^) says that the vector 
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(13) (Cm{^ + 2B-^^7i) : 1 < m < L, fx e K 

is orthogonal to each member of the above orthonormal basis of C°^. Hence, the 
vector in the expression (|13D is zero. In particular, Cm{0 — fo'^ all m, 1 < m < L. 
That is, Cmp = 0, l<m<L, Z'^. Therefore, / = 0. This ends the proof. □ 
The splitting lemma can be used to decompose an arbitrary Hilbert space into 
mutually orthogonal subspaces, as in 0. We will use the following corollary later. 

Corollary 2.3. Let : 1 < I < L, k E Z^} be an orthonormal basis of a sepa- 
rable Hilbert space Ti. Let Hr, < r < a — 1 be as above and satisfy (p!OD. Define 

Flu = t E hl^^^_^,E^p- 0<r<a-l, 1<1<L, keZ''. 

Then {F[^ : 1 < I < L, k e Z'^} is an orthonormal basis for its closed linear span 
W andn = ®'i.zlW. 

Proof: Let yji, (y92; • • • ? V^l be functions in L^(M'^) such that — /c) : 1 < / < L, 

k E Z'^} is an orthonormal system. Let V = 'sp{a^^'^(pi{A ■ ~k) : I, k}. Define a 
linear operator T from the Hilbert space to H by T(a^^'^Lpi(A ■ —k)) = Ei^k- Let 
f[ are as in (^. Then, T{f[{- — k)) = F[f^. Now the corollary follows from the 
splitting lemma. □ 

3. Construction of Multiwavelet Packets 

Let {Vj : j e Z} be an MRA of L^{R'^) of multiplicity L associated with the 
dilation matrix A. Let {(pi : 1 < / < L} be the scaling functions. Since ipi,l < I < L 
are in Vq C Vi and {a^^'^(pj{A ■ —k) :1 < j < L,k E Z"'} forms an orthonormal basis 
of Vi, there exist {hijk : k e Z'^} e f{Z'^) for 1 < /, j < L such that 

L 

'^li.^) = E E hijkO^/'^Lpj^Ax - k). 
Taking Fourier transform, we get 

= tj: h,,,a-'/'e-<^-'^'')0,{B-'O 

(14) = thiAB-'O0AB-'O, 

i=i 

where hij{^) = Yl a^^^^hijke~''^^'''\ and hij is 27rZ'^-periodic and is in L^(T°'). Let 
HoiO be the L x L matrix defined by 

We will call i^o the low-pass filter matrix. 
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Rewritting ([1^ in the vector notation (|]) and (|^), we have 



(15) l>(0 = Ho{b-'oHb-'0- 

Let Wj be the wavelet subspaces, the orthogonal complement of Vj in V^+i: 



Properties (Ml) and (M3) of Definition L2 now imply that 

W,±Wj,, j^f 

and 

(16) feWj^ f{A-^-) e Wo. 

Moreover, by (M2), L^{R'^) can be decomposed into orthogonal direct sums as 

(17) L\R'') = ^Wj 

jez 

(18) = Vo®[^W,). 



By Lemma |2.1| and equation ([T^, we have (for I < l^j < L) 



E {Y.h,^{B-\^ + 2k'K))^^{B-\^ + 2k^))] 



Now, using (H), we have 



L L 



Sji = E E E ^.m(5-ie + 25-Vvr)/izn(S-ie + 25-Vvr) 

^&Kb ra=l n=l 

■ E {^™(fi"'(e + 2/i7r) + 2kTi)^n{B-\i + 2/i7r) + 2A;7r)|, 
where Kb is a set of digits for B. Using Lemma ^?T] again, we get 

(19) 8ji= E E h,^{B-^i + 2B-^fi7r)hUB-'^ + 2fi-V7r). 

H&Kb m=l 

This is equivalent to saying that 

E ^o(e + 2B-^fni)H*{^ + 25-Vvr) = /l for a.e. C 

Equation (|19|) is also equivalent to the orthonormality of the vectors 

(hij{^ + 2B~^i2n) : 1 < j < L, G Kb) , 1 < / < L, ^ G T^. 
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These L orthonormal vectors in the aL-dimensional space C^^ can be completed, 
by Gram-Schmidt orthonormahzation process, to produce an orthonormal basis for 
C^^. Let us denote the new vectors by 

[h\^{i + 2B-^fiTT) : l<j<L,i2e Kb), I < I < L, 1 < r < a - 1, ^ G 

and extend the functions h'^j {1 < r < a — 1, 1 < l,j < L) 27rZ'^-periodically (see 
for the one-dimensional dyadic dilation). 

Denoting hj Hr{^), 1 < r < a — 1 the L x L matrix 

we have 

^ Hr{^ + 2B-^^-K)Hl{i + 25-Vvr) = Srsh for a.e. ^. 
Now, for l<r<a — 1, 1 < I < define 

(20) fI{i)=th\,{B-'i)vAB~'i). 

Since h\j are 27rZ'^-periodic, there exist {h\jf, : k G Z'^} G P{Z'^) such that 

Now, applying the splitting lemma to Vi, we see that {f[{- — k) : < r < a — 1, 
1 < Z < L, /c G Z'^} is an orthonormal basis for Vi. We use the convention 
= fiA <^ < L with hij = h^ij and hi^k = h^ijk- 

The decomposition Vi = Vo® Wq, and the fact that {f^{—k) : I < I < L, k e Z'^} 
is an orthonormal basis of Vq, imply that 

{//■(■ - A;) : 1 < r < a - 1, 1<1<L, k e Z'^} 

is an orthonormal basis for Wq. By ( |16|) and ([T7|) , we see that 

{a^/'^f[{A^ ■ -A;) : 1 < r < a - 1, 1 < / < L, j G Z, A; G Z"^} 

is an orthonormal basis for L^(]R'^). This basis is called the multiwavelet basis and 
the functions {f[ : I < r < a — 1, I < I < L} are the multiwavelets associated with 
the MRA {Vj : j G Z} of multiplicity L. 

For < r < a — 1, by denoting Fr{x) = (/[(a;), /2 (a;), . . . , /£(a;))* and -Fr(0 = 
(/[(O, /2 (0, • • • , we can write (pD and (10) as 

(21) FriO = H,{B-^OHb-'0, < r < a - 1. 

This equation is known as the scaling relation satisfied by the scaling functions 
(r = 0) and the multiwavelets (1 < t < a — 1). 

As we observed, applying splitting lemma to the space Vi = sp{a^^'^(pi{A ■ —k) : 
1 < / < L, G Z"^}, we get the functions fl, < r < a - 1, 1 < / < L. Now, for 
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any n G Nq = NU {0}, we define 1 < I < L recursively as follows. Suppose that 
//", r G No, I ^ I < L are defined already. Then define 

(22) /;+-(a;) = E E Mj^/'fJiAx - k); < s < a - 1, 1 < I < L. 
Taking Fourier transform 

(23) if mo = thUB-'oifjriB-'o- 



In vector notation, (123| ) can be written as 

(24) {Fs+arnO = H,iB-^OFr{B-'0- 

Note that (^) defines /" for every non-negative integer n and every / such that 
1 < I < L. Observe that fi = ^lA < I < L are the scaling functions and 
/;'",l<r<a — 1,1</<L are the multiwavelets. So this definition is consistent 



with the scaling relation (|2T|) satisfied by the scaling functions and the multiwavelets. 

Definition 3.1. The functions {/" : n > 0, 1 <l < L} as defined above will be 
called the basic multiwavelet packets corresponding to the MRA {V, : j G Z} of 
L'^{W^) of multiplicity L associated with the dilation A. 

The Fourier transforms of the multiwavelet packets 

Our aim is to find an expression for the Fourier transform of the basic multiwavelet 
packets in terms of the Fourier transform of the scaling functions. For an integer 
n > 1, we consider the unique "a-adic expansion" (i.e., expansion in the base a): 

(25) n = yUi + /i2a + n^a^ + ■ ■ ■ + Hja^~^, 

where < /ij < a — 1 for alH = 1, 2, . . . , j and fij ^ 0. 

If n can be expressed as in (p5[) then we will say n has a-adic length j. We claim 
that if n has length j and has expansion (^), then 

(26) F„(0 = H^,{B~^OH,AB-\) ■ ■ ■ H,^{B-^O^B-^0 

so that (/")^(0 is the l-th component of the column vector in the right hand side 
of (pBI). We will prove the claim by induction. 

From ( pH ) we see that the claim is true for all n of length 1. Assume it for length 
j. Then an integer m of a-adic length j + 1 is of the form m = ^ + an, where 
< fi < a — 1 and n has length j. Suppose n has the expansion (P3|). Then from 
(P^) and (P^D, we have 



= H^{B-'OF„iB-'0 

= H,{B''i)H,, {B-'O ■ ■ ■ H,^ {B-^^^'^OHb-^'^'^0- 

Since m = fi + an = fi + fiia + fi2a^ H — ■ + fija^ , -^m(0 the desired form. Hence, 
the induction is complete. 



12 BISWARANJAN BEHERA 

The first theorem regarding the multiwavelet packets is the following. 

Theorem 3.2. Let {/" : n > 0, 1 < I < L} be the basic multiwavelet packets con- 
structed above. Then 

(i) {/"(■ — k) : < n < a^^^ — 1, 1 < I < L, A; G Z'^} is an orthonormal basis 
ofW,,j>^. 

(ii) {/"(■ — /c) : < n < a-' — 1, 1 < / < L, fee U^} is an orthonormal basis of 

y,. J > 0. 

(iii) {/"(• - A;) : n > 0, I < I < L, k E U^} is an orthonormal basis of L'^iW^). 

Proof: Since {/" :l<'n,<a — 1, 1 <l < L} are the multiwavelets, their U^- 
translates form an orthonormal basis for Wq. So (i) is verified for j = 0. Assume 
for j. We will prove for j + 1. By assumption, the functions {fP{- — k) : a^ < n < 
a^^^ — 1, 1 < I < L, k E li'^} is an orthonormal basis of Wj. Since / G Wj <^ 
f{A-) G Wj^i, the system of functions 

{a^/^fPiA ■ -k) ■.a^<n< a^+^ - 1, 1 < / < L, keZ'^} 

is an orthonormal basis of Wj+i. Let 

En = sp{a^/^fJ\A --k): 1<1<L, ke Z'^}. 

Hence, 

(27) W,^, = En. 
Applying the splitting lemma to we get the functions 

(28) g^'\x) = E E h\„,,a'" fli^^ ' k) ( < r < a - 1, 1 < / < L) 

"1=1 fcgZ'* 

so that {g]^'^{- — k) : < r < a — 1, 1 < I < L, A;G Z'^} is an orthonormal basis 
of En. But by (p2|), we have 

yi — Jl 



This fact, together with (|27|) , shows that 

Ijr+anj^. _ ^) . Q < r < a - 1, 1 < / < L, G Z^, o^' < n < a-'+^ - 1} 
= {/r(- - k) : a^^^ <n< a^+^ - 1, 1 < / < L, A; G Z"'} 

is an orthonormal basis of W^j+i- So (i) is proved. 

Item (ii) follows from the observation that Vj = Vq® Wq ® ■ ■ ■ © W^j-i and (iii) 
follows from the fact that vWj = ^^(M"'). □ 
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4. Construction of Orthonormal Bases From the Multiwavelet 

Packets 

We now take all dilations by the matrix A and all Z'^-translations of the basic 
multiwavelet packet functions. 

Definition 4.1. Let {fp : n > 0, 1 < Z < L} be the basic multiwavelet packets. 
The collection of functions 

r = {a^/Vr(^^ ■-k):n>0, 1 < I < L, j e Z, k e Z'^} 

will be called the general multiwavelet packets associated with the MRA {Vj} 
of L'^{R'^) of multiplicity L. 

Remark 4.2. Obviously the collection V is overcomplete in L^(R'^). For example 

(i) The subcoUection with j = 0, n > 0, 1 < I < L, k E Z'^ gives us the basic 
multiwavelet packet basis constructed in the previous section. 

(ii) The subcoUection with n = 1, 2, . . . , a — 1; 1 < I < L, j e k e Z'^ is the 
usual multiwavelet basis. 

So it will be interesting to find out other subcollcctions of V which form orthonor- 
mal bases for L^(M'^). For n > and j e Z, define the subspaces 

(29) = sp{a^/Vr(^^ --k)- 1<1<L, keZ'^}. 

Observe that 

a-1 



U^ = Vj and ^Uj=Wj, j e Z. 



r=l 



Hence, the orthogonal decomposition Vj+i — Vj® Wj can be written as 



a-1 



t^.i = e^;- 



r=0 



We can generalize this fact to other values of n. 
Proposition 4.3. For n> and j e Z, we have 

a-1 

(30) [/jYi = 0^r^'^- 

r=0 

Proof: By definition 

C/jYi = sp{a'^fP{A^+^ --k): 1<1<L, ke Z'^}. 

Let 



Ei^k{x) = a'^fPiA^^^ ■ -k), for 1 < / < L, ke Z^. 
Then {Ei^k ■ I < I < L, k e Z'^} is an orthonormal basis of the Hilbert space t^JYi- 
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For < r < a - 1, let 



m=i/3eZ'* 



and 



rc = sp{F[k ■■ 1<1<L, ke Z^}. 

Then, by Corollary |2.3| we have 



a-l 



r=0 



Now 



E E him,l3-AkFm,,p{x) 
L 

E E h\,m,aFm,Ak+a{x) 



Therefore, 



and 



a'E E 



-^r jjan+r 



^]Vi = ©t^; 



an+r 



r=0 

□ 

Using Proposition ^]3| we can get various decompositions of the wavelet subspaces 
Wj, j > 0, which in turn will give rise to various orthonormal bases of L^(M'^). 

Theorem 4.4. Let j > 0. Then, we have 

a-l 



r=l 
a2-l 

= 0^-1 



(31) 
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K, 



l<3 



r=aJ> 



where UJ is defined in (pOf). 

Proof: Since Wj = Q^lj f/J, we can apply Proposition |4.3| repeatedly to get (|3lD . 

□ 

Theorem |4.4| can be used to construct many orthonormal bases of L^(M'^). We 
have the following orthogonal decomposition (see (|l3)): 



Vo®Wo®Wi® W2 



For each j > 0, we can choose any of the decompositions of Wj described in (|3TD . 
For example, if we do not want to decompose any Wj, then we have the usual 
multiwavelet decomposition. On the other hand, if we prefer the last decomposition 
in (|31|) for each Wj, then we get the multiwavelet packet decomposition. There are 
other decompositions as well. Observe that in (|31|), the lower index of f/"'s are 
decreased by 1 in each succesive step. If we keep some of these spaces fixed and 
choose to decompose others by using (^), then we get decompositions of Wj which 
do not appear in (PT|). So there is certain interplay between the indices n G Nq and 

jez. 

Let S* be a subset of Nq x Z, where Nq = N U {0}. Our aim is to characterize 
those S for which the collection 

Vs = {a'^fnA' --k): 1 < I < L, keZ^ (n, j) G s] 

will be an orthonormal basis of L'^{W^). In other words, we want to find out those 
subsets 5 of No X Z for which 



(32) 



f/; = L2(R'^). 



By using (^) repeatedly, we have 



a-l 



(33) 



ur.r 

r=0 

a(n+l)-l 



a(n+l)-l 



a-l 

0K 

.s=0 



ar+s 

i-2 







i-2 







jjr 
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Let In,j = {r G No : a^n < r < a^(n + 1) - 1}. Hence, 
(34) = U^o- 

That is, 

{n,j)eS {n,j)eS r-6-f(nj) 

But we have already proved in Theorem p.2| that 

rgNo 



Thus, for (|3^) to be true, it is necessary and sufficient that {Inj ■ ('^jj) G 5} is a 
partition of Nq. We say {Ai : / G /} is a partition of Nq if Ai C Nq, Ai's are pairwise 
disjoint, and Ui(zjAi = Nq. We summarize the above discussion in the following 
theorem. 

Theorem 4.5. Let {/" : n > 0, 1 < I < L} be the basic multiwavelet packets and 
5 C No X Z. Then the collection of functions 

[aifJ'^A^ --k): 1 < I < L, keZ^ (n,j) G 

is an orthonormal basis of L'^{M.'^) if and only if {I^j '■ {n,j) & S} is a partition of 
No. 

5. Wavelet Frame Packets 

Let Hhe a separable Hilbert space. A sequence {xk : k E Z,} oi H is said to be a 
frame for Ti. if there exist constants Ci and C2, < Ci < C2 < cxd such that for all 

(35) C-illxf < J2\{x,^k)\^ <C2\\xf. 



The largest Ci and the smallest C2 for which (^3p holds are called the frame 
bounds. 

Suppose that <l> = . . . ,(p^} C L'^{R'^) such that Wi- - k) : 1 < I < N, 

k G Z°'} is a frame for its closed linear span «S'(<I>). Let ip'^ , ip"^ , . . . ^ip^ be elements 
in 5'($) so that each ip^ is a linear combination of </?'(■ — /c); 1 < / < L, G T/. A 
natural question to ask is the following: when can we say that {ip^{- — k) '■ 1 < j < N, 
k G Z'^} is also a frame for 5'($)? 

If G S{^), then there exists {pjik : A; G Z'^} in such that 

TV 

^j(x) = XI E Pjik<^\x - k). 

In terms of Fourier transform 

N 

1=1 
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N 

1/ 



(36) = EPjim'iO {1<J<N) 



1=1 

where Pji{^) = ^ Pjike~'^'''^l Let P(0 be the N x N matrix: 

V / l<j,l<N 

Let 5* and T be two positive definite matrices of order A^. We say S* < T if 
{x,Sx) < {x,Tx) for all x G M^. The following lemma is the generalization of 
Lemma 3.1 in 

Lemma 5.1. Let ip\%l)^ for 1 < I < N , and P{C,) be as above. Suppose that there 
exist constants Ci and C2, < Ci < C2 < oo such that 

(37) C^I < P*iOP{0 < C2I for a.e. ^ e 

Then, for all f G ^^(R"'), we have 

(38) 

N o N „ TV „ 

^lE E|(/,^'(--fc)>| < E EK/,^'(--fc)>| <C2E EK/,^'(--fc)>| • 

l=lk&'i 1=1 kel'^ l=lk&'i 

Let A be a dilation matrix, B = A* and a = \ det A\ = \ det B\. Let 

(39) Ka = {ao,ai, . . . ,aa-i} 
and 

(40) KB = {(3o,f3,,...,(3a-i} 

be fixed sets of digits for A and B respectively. For < r, s < a — 1 and 1 < I, j < L, 
define for a.e. ^ 

(41) SlfiO = (5z,a-5e-K?+2S-'/5^-'"'-) . 
Let 

(42) E'^iO={£IfiO)^^ 
and 

(43) E{0 = (E^'iO) 

So E{^) is block matrix with a blocks in each row and each column, and each 
block is a square matrix of order L, so that E{C,) is a square matrix of order aL. We 
have the following lemma which will be useful for the splitting trick for frames. 

Lemma 5.2. (i) If u e Ka, then ^ ^-i2n{B-^,.,u) ^ 
(ii) The matrix E{C,), defined in (^), is unitary. 



1<IJ<L 
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Proof: Item (i) is the orthogonal relation for the characters of the finite group 
Tf' jBTj^ (see Observe that the mapping 

is a character of the (finite) coset group /B7/. If = (i.e., if z/ G AT/), then 
there is nothing to prove. Suppose that 7^ 0, then there exists a yu' G Kb such 

that e"*^'^'^^ ^ 7^ 1. Since Kb is a set of digits for B, so is Kb — /u'. Hence, 
Now 

^ ^-i2n{B-^fi,u) ^ g-j27r(B-V','^) . ^ g-j27r(B-l(AJ-/x'),i.) 

= g-i2^(i?-V',^) . ^ g-i27r(i?-V,i')^ by (H). 

Therefore, 

E e-'2.(B-V,-) = 0, since e-''"" ^' ^ 1. 
AteA's 

To prove (ii), observe that the (r, s)-th block of the matrix E{^)E*{^) is 

t=o 

The (/, j)-th entry in this block is 

t=0 m=l 

a-1 L 



=Om=l 

E 5;m5,,nEa-'e-^(«+2^"^-'"^-"=) 



m=l t=0 
L 

= E SimSjmSrs, (by (1) of the lemma) 

m=l 
= SijSrs- 

This proves that E{^)E*{^) = I. Similarly, E*{^)E{^) = I. Therefore, E{^) is a 
unitary matrix. □ 

6. Splitting Lemma for Frame Packets 

Let {^i : I < I < L} he functions in /.^(R'^) such that {</)/(■ - k) : 1 < I < L, 
k G Z'^} is a frame for its closed linear span V. For < r < a — 1 and 1 < I < L, 
suppose that there exist sequences {h^jf, : k G Z''} G t^il/). Define fl as in (^) and 
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That is, 

(45) f[{x)=t E h^i,,a'/'vAAx-k). 

Let Hr{^) be the matrix defined in (^. Let and Kb be respectively fixed sets 
of digits for A and B as in (^) and (^D|). Let H{^) be the matrix 

(46) H{0 = (HrU + '^B-'Psn)) 

V / 0<r,s<a-l 

H{C,) is a block matrix with a blocks in each row and each column, and each block 
is of order L so that H{^) is a square matrix of order aL. 

Assume that there exist constants Ci and C2, < Ci < C2 < oo such that 

(47) CJ < H*{^)H{^) < C2I for a.e. ^ G T'^. 
We can write f[ as 

= E E E (^^ - - ^fc), by (D 

j=is=Ofcezrf 

L a—1 

= E E E ^[,j,a.+Afc</'i (a; - k), 

i=ls=OA:eZd 

where 

(48) (pf{x) = a^/^(fj{Ax - as), < s < a - 1. 
Taking Fourier transform, we obtain 

= EEW(0(v^f)'(0, 

3=1 s=0 

where p[/(0 = E ^F,i,a,+Afce"^<^'''^ Define 

(49) 

and 



0<r,s<a-l 



(50) P(0 = (P''^(0 
Claim: 

(51) H{i) = P(i?0^(0 
where E{^) is defined in (|l])-(|3D. 
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Proof of the claim: The (r, s)-th block of the matrix P[BC,)E{C,) is the matrix 



t=o 

The (/, j)-th entry in this block is equal to 

t=Om=l 
a-1 L 



i(BC,A:> -1/2 -i(C+2B-i/3«7r,Qt) 



t=Om=l 
a-1 

E E Kj,o^t+Ak(^ "'"'"^'^ 



Now, the (Z, j)-th entry in the (r, s)-th block of H{C,) is 



-1/2 



a-1 



,-1/2 



t=0 

a-1 



E E K,j,at+Ak^ 



t=0 



So the claim is proved. In particular, we have 



(52) 



Since E{^) is unitary by Lemma |]2], and P*{B^)P{B^) are similar 

matrices. Let A(^) and A(^) respectively be the minimal and maximal eigenvalues 
of the positive definite matrix H*{C,)H{^), and let A = inf A(^) and A = sup A(^). 

(It is clear from (|5lD that A(^) and A(^) are 27rZ'^-periodic functions.) 

Suppose < A < A < oo. Then we have, by (|^) (in the sense of positive definite 
matrices) , 

XI < H*{^)H{^) < AI for a.e. ^ G T'^ 

which is equivalent to 

XI < P*{^)P{^) < AI for a.e. ^ G T'^. 
Then by Lemma f).l[ for all g G L^(M"'), we have 



a-1 L 

AEE E 



is). 



k))' < EE E \{9jn--m 

s=Ol=ll 



(53) 



a-1 L 



< AEEE {9,vl'\--k) 
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Since 



(s) 

where (pl is defined in 



L 2 a— I L 



(54) E E \{9,a'/'MA--k))\'=EE E (9,vr\-k) 

which follows from (PSI), inequality (El) can be written as 



\{9,a'/'MA--k))\' < J:tj:\{9,fn--k))f 

(55) < AE E \{9,a'/'^iiA.-k))\\ 

This is the splitting trick for frames: the 74~^Z°'-translates of the L dilated func- 
tions ipi{A-), 1 < I < L, are 'decomposed' into Z'^-translates of the aL functions , 
0<s<a-l, 1<1<L. 

We now apply the splitting trick to the functions {/f : I < I < L} for each s, 
0<s<a — Ito obtain 

^tE\{9,a'^'fnA--k))\' < EEE \{9jri-k))f 
l=ikeZ'i r=oi=ikeZ'i 

(56) < AE E \{9,a'/'fnA--k))\\ 

1=1 k&d 

where < r < a — 1 are defined as in (^Sj) (/f now replaces ipi): 

(57) f^ix) = E E htj,a'/'fJ{Ax - k); < s < a - 1, 1 < I < L. 

j=i fcez<* 

Summing (0) over 0<s<a — 1, we have 

AEE E K^7,«^/Vf(^--A:)>r < EEE E l(^7,/r(--fc))l' 

s=Ol=lkeZ'^ s=0r=0l= 



Using (|55|), we obtain 



< AEE E \{9,ay^ft{A.-k))\ 

S = Ol = lke:Z^ 



X'tj: \{9,a'/'MA'--k))\' < EEEE \{9,fri--k))f 

l=lk€Z'i s=0 r=0 1=1 keZ"^ 

(58) < A^E E \{9,a'/'^iiA'--k))\\ 

1=1 k€Z^ 



n 



Now as in the case of orthonormal wavelet packets, we can define for each 
> and 1 < / < L (see (||) and (|26|)). In order to ensure that are in L'^{R'^), 
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it is sufficient to assume that all the entries in the matrix H{C,), defined in (^61), are 
bounded functions. Comparing { ^7\j and (p2|), we see that 

{fr : < r, s < a - 1} = {/;+'^'^ : < r, s < a - 1} = {fl" : < n < - 1}. 

So ( pSD can be written as 



L 



< A'E E \{9,a'/'MA'--k))\\ 



1=1 



By induction, we get for each j > 1 



A^'E E K^?,«^'/V/(^^'--A:)>r < e'e E //"(■- 

(59) < A^'E E K^7,a^'/VK^^'--A:)>| 



We summarize the above discussion in the following theorem. 

Note: {/" : n > 0, 1 < I < L} will be called the wavelet frame packets. 

Theorem 6.1. Let {^i : I < I < L} C L'^{R'^) be such that - k) : 1 <l < L, 
k G Z'^} is a frame for its closed linear span Vq, with frame bounds Ci and C2. 
Let H{^), Hr{^), X and A be as above. Assume that all entries of Hj.{^ + 2B~^(3s7r) 
are bounded measurable functions such that < A < A < cxd. Let {/" : n > 0, 
I < I < L} be the wavelet frame packets and let Vj = {f : f{A~^-) G Vq}. Then for 
all j > 0, the system of functions 

{fP{- - A;) : < n < a-'' - 1, 1 < I < L, A; G Z^} 

is a frame ofVj with frame bounds X^Ci and A!' €2- 

Proof: Since — k) : 1 < / < L, A; G Z"'} is a frame of Vq with frame bounds 
C\ and C2 , it is clear that for all j 

{a^''^^i{A^ --k): 1 <l < L, k eZ'^} 

is a frame of Vj with the same bounds. So from (0), we have 



(60) A^Cill^^f < e'e E \{9jn--k))f <A'C2\\9r for all ^? G IS-- 

□ 

In Theorem ^]2| we proved that the basic multiwavelet packets form an orthonor- 
mal basis for L^(M°') = Ul^-. An analogous result holds for the wavelet frame packets 
if the matrix H{^), defined in (|46D, is unitary. 
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Before proving this result let us observe how the space ^j>oVj looks like. Let 
Vo = spWii- -k):l<l<L,ke Z^}, = {/ : f{A~^-) E Vo} and Vj C V,+i. 
Let W = UVj. Then it is easy to check that f e W ^ f{- - A~^k) G W for all 
j E 1j and k G Z'^. We claim that elements of the form A~^k are dense in R''. For 
K = {ki, k2, . . . ka} a set of digits for A, define the set 

Q = Q{A, K) = l^x eR'^ : X = J2 ^'^^v ^ 

In the above representation of x, e^-'s need not be distinct. We have 

\\A-^x\\ < Ca^\x\\, X G 

where C is a constant and < a < 1 (see [|r^, Chapter 5]). Therefore, the series 
that defines x is convergent. 

For X = {xi,X2, . . . , Xd) G M'^, \\x\\ = + |x2p + ■ ■ ■ + |a;rfp)2 . 

The set Q satisfies the following properties (see [p!0| ): 

(i) Q = UUA-HQ + k^) 

(ii) Uk&iQ + A;) = 

(iii) Q is compact. 

Let e > and y E Q. We first show that there exist j G Z and /c G Z"^ such that 
\\y — A^^kW < e. From (i) we have 

a 

Q = [jA-\Q + h) 

i=l 
a 

= u-^- 

i=l 

a a 

= [j[j{A-'Q + A-^km + A-'ki). 

i=l 771=1 

Hence, for any j > I and any y & Q, there exist yj G Q and li,l2, ■ ■ ■ Jj ^ K such 
that 

y = A-^yj + A-Hj + A-^^-^^lj^i + ■■■ + A-Hi. 

Therefore, 

\\y-A-^m + Al,^, + --- + A^'\}\\ = \\A-^y,\\ 

< Ca^y^W 

< C'a^ (as Q is compact) 

< e, choosing j suitably. 

Now if y G M*^, then by (ii) y = + p ioi some y^ E Q and p G U^. For yQ G Q, 
there exist j > and /c G Z*^ such that ||?/o ~ A^-'A;]! < e. That is, 

\\yQ+p- A-^{k + A^p)\\ < e 
^ \\y-A-^{k + A^p)\\<e. 



a 

U A-\Q + k^)+h 

m=l 
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So the claim is proved. 

We have proved that W is invariant under translations by A~^k and these elements 
are dense in W^. Therefore, FT is a closed translation invariant subspace of L^(M'^). 
Hence, W = I/|(]R'^) for some E cR'^ (see where 

Ll{R^) = {f e L\R'^) : supp f C E}. 

Now let 

L 

Eo = [j\jB^{suppifi). 

1=1 j>0 

Claim: E = Eq a.e. 

To prove the claim we will follow Theorem 4.3. Since (fi{A^-) E Vj G W, the 

function = jj0i{B~^-) eW = {f : f E W}. Therefore, B^{supp ^i) = 

supp C E for all j > and 1 < I < L, which implies that Eq C E. 

Let Ei = E\Eo.We have 

(61) fEV,^f= tmiB-^O^iiB-'O, 

1=1 

for some 27rZ''-periodic functions mi E L^(T'^). 

Hence, (0) implies that f = on Ei for all f E Vj and hence, for all / E UVj = W. 
Taking closure, we obtain that / = on ii^i for all / E W . But W is the set of all 
functions whose Fourier transform is supported in E. Since Ei C E, we get that 
£^1 = a.e. Therefore, E = Eq a.e. □ 

Theorem 6.2. Let {ipi{- - k) : 1 < I < L, k E Z"^} C L'^{R'^) be a frame for its 
closed linear span Vq, with frame bounds Ci and C2 and let Vq C Vi, where Vj = 
{/ : /(A^-'-) E Vo}. Assume that H{^) is unitary for a.e. C,- Then {/"(■ — k) : n > 
0, 1 < / < k E Z'^} is a frame for the space ^j>oVj with the same frame bounds. 

More generally, let S = {(n,j) G No x Z} be such that U(nj)e5-^".i ^■^ partition 
o/Nq. Then the collection of functions {a^^'^fl^{A^ ■ —k) : 1 < I < L,(n,j) E S, 
k E Z*^} is a frame for ^j>oVj with the same bounds Ci and C2. 

Proof: Since H{C,) is unitary, A = A = 1 so that the inequalities in (|59D are 
equalities, and from ( pOD we have 

(62) Ci||^7f < e'EE \{9Jri-~m'<C2\\9f ioT^WgEVj. 

n=Ql=lk&d 

Now let h E Uj>o^'. Then there exists hj E Vj such that hj —^hasj^ 00. Fix 
J, then for j < j', we have from ( p^ ) 

e'e e \{hj',fn--k))f<c2\\hA'- 

n=Q l=lkeZd 

Letting j' 00 first and then j — > 00, we have for all h E ^j>oVj 
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(63) EE E \{hjn--k))\'<C2 

n>Ol=lkeZ'i 

To get the reverse inequality we again use (|6^): 

Ci\\h,r < e'e e \{h,jn--k))\' 



n=0 l=lk&'^ 

e'e e \{h,-hj-{.-k)) + {hjr{--k))t 

n=0 1= 



Therefore, 

n=0 1=1 keZ'i ' 

+( EE E \{hjn--k))\y 
< ci"\\h,-h\\ + r^Y:j: E \{hjn-~m")\ ^^jm 



n=0 [=\k&L'i 

Taking j — > cxo, we get 



c.\\hf< \{h, fn- - k))f 

n>Ol=lk£Zd 



as 



for all h G UV^-. So the first part is proved. 

Now let = sp{a^/^fp{A^ ■ -k) : 1 < I < L, ke Z'^}. Then we can prove 
in the orthogonal case (see (|3^)) that 

where ^ is just a direct sum not necessarily orthogonal, and = {r G No : a^n < 
r < {n + 1) — 1}. Now, since H{^) is unitary, we have A = A = 1 and hence (^) 
is an equality. Therefore, 

E E K^,a^/Vr(^--fc)>r= EE E \{9jr^''{--k))\\ 

1=1 k&d r=Ol=lkeZ<i 

From this we get 

EE K^,«'/Vr(^'--^)>r=EEE E (9Ji^^''^'^^\--k) 

a2(n+l)-l L 

= E EE \{9J[i--k))f. 

r=a?n 1=1 k&Z'^ 
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Similarly, 

L „ ai{n+l)-l L 

E E \{g,a=''fr{A^--k))\'= E EE \{9J[i--k))\ 

(64) = E E E \{9j[{--k))\\ 



2 



From the first part of the theorem, we have for all / G UVj 

Ci\\fr< EE E \{fjn--k))\'<c. 

But, the set S is such that IJ(nj)6S-^rij = Nq. Therefore, 

Cill/f < E E E E \{fj[{--k))f<c. 

Using (p^), we get 

Ci||/f< E EE \{f,a^^'fnA^--k))\'<C2 

{n,j)esi=ikeZ'i 



for all / G UVj-. This completes the proof of the theorem. □ 
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